Generating sub-TeV quasi-monoenergetic proton beam by an ultra-relativistically 

intense laser in the snowplow regime 
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Snowplow ion acceleration is presented, using an ultra-relativistically intense laser pulse irradi- 
ating on a combination target, where the relativistic proton beam generated by radiation pressure 
acceleration can be trapped and accelerated by the laser plasma wakefield. The theory suggests that 
sub-TeV quasi-monoenergetic proton bunches can be generated by a centimeter-scale laser wakefield 
accelerator, driven by a circularly polarized (CP) laser pulse with the peak intensity of lO^^W/cm^ 
and duration of 116fs. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw, 52.59.-f 



Laser-driven ion acceleration has drawn attention for 
many applications, e.g., proton cancer therapy fast ig- 
nition of thermonuclear fusion by protons [2| , conversion 
of radioactive waste J3| , high energy physics accelerator 
Qand astrophysics All these applications may be 
enabled by the introduction of near-future high-intensity 
lasers (i.e. ELI), which will be capable of producing 
pulses with intensity 10^^ — 10^'' W/cm^ (6|. 

Radiation Pressure Acceleration (RPA)[3l accelerates 
ions efficiently and theoretical studies show that GeV 
proton beam may be generated by an ultra-intense laser 
with intensities above lO^^W/cm^ [1, [13, Eli. Ahhough 
the acceleration field is over tens TeV/m, the accel- 
eration length is normally shorter than hundreds mi- 
crons. Because of the energy scaling, it is difficult to 
further increase the ion to very high energies (i.e. TeV ). 
Quasi-monoenergetic GeV electron beam can be gener- 
ated by a centimeter-scale laser plasma wakefield acceler- 
ator(LPWA) 12|; however, it is not easy to capture and 
continuously accelerate the nonrelativistic ions. Shen et 
al. [isj reports a laser pulse with an ultrarelativistic inten- 
sity can excite an electrostatic field, which can capture 
protons from underdense plasma and accelerate them to 
tens GeV. Recently Yu et al. 14| found that the under- 



dense gas behind a thin foil can accelerate the proton 
beam to 60GeV by driving high-amplitude electrostatic 
fields moving at a high speed. There is an important 
question whether it is possible to generate TeV-level pro- 
ton beams by the laser plasma wakefield accelerator. 

This letter reports an ultra-relativistically intense laser 
can excite the plasma wakefield that traps the relativis- 
tic proton beam and accelerates it over a few centime- 
ters. An analytic model is derived for proton acceler- 
ation that extends the applicability of Esarey's classic 
plasma wakefield theory jl5| and has been confirmed by 
our PIC simulations. Our scaling law shows that sub- TeV 



quasi-monoenergetic proton bunches can be generated by 
a circularly polarized (CP) laser pulse with the intensity 
of 10^'^ W/cm^ and duration of 116 fs. A combined tar- 
get in FigH] (a) shows a planar hydrogen target that is 
followed by underdense heavy-ion background gas with a 
charge-to-mass ratio of 1/3. When the foil thickness D is 
equal to j^^o-o^h a mono-energetic proton beam is gen- 
erated from a moving double layer Figl^Ja) suggests 
that protons and electron layer are very close. The rel- 
ativistic protons can quickly overrun the laser front that 
results in shorter dephasing length and lower proton en- 
ergy. If a thinner foil satisfies 
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the ponderomotive force quickly pushes electrons to the 
rear of the foil (see Figl^Jb)), the electrons that are piled 
up there reaches density far greater than jUc. Mean- 
while protons in the foil are accelerated by the separated 
charge field. Here lo,ne,nc, and A/ are the plasma skin 
depth, plasma density,critical plasma density, and laser 
wavelength, respectively, oq = eA/meCOic is the normal- 
ized laser amplitude and A is the laser vector potential. 
7 = (l+/Af/1.37xl0i^)i/2, /is laser intensity. A; in units 
of micron. The electron-layer in rear of the foil is pushed 
out of the foil by the laser pulse soon before double-layer 
(consists of electrons and protons) is formed. The layer 
runs in the lower density gas and is further pushed by the 
laser pulse like light-sail to generate " Snow-plow" effects 
on enhancement of electron number. The electrons in gas 
plasma drew away from the foil generate a large electric 
field, then the protons that have been accelerated to GeV 
level start long march in gas by the wakefield accelera- 
tion. This acceleration scheme is named the snowplow 
regime. We also notice a connection of the above plasma 
acceleration with "snowplow" acceleration in Ref. [ig} . 
In order to grasp the main physics and to construct a 
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FIG. 1: (color online) (a) Acceleration scheme, the initial den- 
sity of the hydrogen foil no/ric = 20, thickness D = 0.5A;; the 
plasma density of the tenuous gas Ue/ric = 0.01, the length of 
the gas plasma is 12000Ai; (b)Sketch map for the snowplow 
process, the dynamical density of ions is rii and electron den- 
sity in the snowplow layer is ricom,- The position B indicates 
the laser pulse front and A is an arbitrary point for calculating 
the local electrostatic field in the snowplow region. 



physical model , in this letter one-dimensional(lD) the- 
ory is derived and testified by ID and 2D PIC simula- 
tions. 

Considering the classic ID plasma wakefield driven by 
a laser pulse two typical features should be no- 

ticed in the ultra-relativistic intensity : (i) Lengthening 
of the plasma wavelength Is — (7nc/rie)^/^A;; 
(ii)Enhancement of the nonlinear wavebreaking field 
Eb — 37r(7ne/nc)^/^i?o, where the critical plasma den- 
sity He = TTTTieC^/e^A^. The highly nonlinear wake- 
field wavelength Is and wavebreaking field Eb are in- 
creased by a factor of 7^/^ over the non-relativistic ver- 
sion, where the electric field strength is normalized by 
Eq = rrieUJic/e. In this regime we adopt a CP laser 
pulse with a wavelength A; = and a peak intensity 
Iq = 1.7125 X lO'^^W/cm^ , corresponding to peak dimen- 
sionless laser amplitude oq = 250. When the foil thick- 
ness D satisfies Eq.(l), the relativistic protons can be 
captured by the excited wakefield in the underdense gas 
and be stably accelerated in the snowplow regime over a 
long distance. An analytic model is developed to estimate 
the acceleration field, dephasing length, pump depletion 

length, and maximum proton energy. The acceleration 

3/2 

length in snowplow regime is proportional to Oq and 
maximum proton energy scales with Uq. It shows that sub 
TeV monoenergetic proton beam is generated in the laser 
intensity of 10^'^ W/cm^. To our best knowledge, it is the 
best efficient proton acceleration mechanism in terms of 
acceleration length and maximum energy. 

First we carried out simulation using a fully relativis- 
tic particle- in-cell code (KLAP) [1, The laser pulse 
has a trapezoidal shape longitudinally with 15 A; flat 
top and 20A; ramp on the rise side. The ramp profile 
is a = aosin^^nt/ 40). At time t=0 it is normally inci- 
dent from the left on a uniform, fully ionized hydrogen 
foil of the thickness D = 0.5/xm and normalized density 



N = no/ric = 20. Behind the hydrogen foil is the tenu- 
ous heavy-ion-gas with a charge-to- mass ratio of 1 /3 and 
a plasma density of n,, = 0.0 Inc. Here the simulation box 
has a size of 12000 A; in x plane and a space resolution 
of 20 cells/ A;. Each cell is filled with 10 macroparticles 
for gas plasma and 20000 macroparticles for foil plasma. 
The initial electron and ion temperature are O.leV. 

After the proton beam is accelerated to GeV level in 
the laser foil interactions, the laser pulse passes through 
the foil and piles electrons in the underdense gas plasma 
as FiglT](b) shows, which is one of most important and re- 
markable features in the snowplow regime. Here the elec- 
tron layer and the positive electrostatic field are called 
the snowplow layer and snowplow field, respectively. The 
snowplow layer can continuously trap electrons in the 
snowplow process, so the acceleration field exists stably 
for a long distance. In the snowplow region, it is found 
that the length of the snowplow region is equal to Is = 
('-fnc/rieY^'^Xi and the snowplow velocity approximately 
equals to the group velocity of laser pulse in the under- 
dense plasma. First we estimate the electrostatic field 
Ea at an arbitrary point xa in the snowplow region. The 
component from ions is E^^ = AireniXA — 4ireni{ls — xa)- 
The contribution from the snowplow layer at position xa 
is Ej^ — ±ll^qsl^. The longitudinal electrostatic field at 

XA is Ea = E'^ + E'^ = Anenei^XA — ^s/2). Therefore, 
the maximum snowplow field at the laser front xb is 

EB^M7ne/n,)^/^Eo. (2) 

Since the injected GeV protons are still slower than the 
laser pulse at the beginning, the ion beam should be in- 
jected into the front side of snowplow region. Afterwards 
protons are gradually accelerated and overrun the laser 
pulse ( see Fig|4lja)). The optimized distance between 
injected proton beam and laser pulse front is defined by 
hnj- It is approximately jls in our simulation. Once 
the proton beam is injected into the wakefield, its maxi- 
mum energy depends on the laser pump depletion length 
and proton dephasing lerigth. The laser pump depletion 
length Lpd is estimated [l5j: E^ ■ Lpd — Ef ■ Li, where 
Ex is the longitudinal electrostatic field driven by laser 
pulse in the snowplow wakefield and Li is the length of 
laser pulse. Assuming the group velocity of laser pulse 
is extremely close to light velocity in free space, pump 
depletion length can be defined as 
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(3) 



where Vetch is the depletion velocity of laser pulse. The 
etching velocity of laser pulse reads Vetch — aon° 
shows the etching velocity of laser pulse scales as n,,/ a^. 

The dephasing length of proton is given approximately 
by Ldp ~ ■ Vp, where kn, is the injected 

length between proton beam and laser pulse front, Vg 
is the laser group velocity, Vp is the proton velocity. It 
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FIG. 2: (Color online) Electron density ne/nc(blue solid 
line), proton density np/nc{the red solid line) and longitudinal 
electrostatic field eE^/meUJic (the black dashed line).(a)When 
D = n^'^O'^' double-layer is formed at t=50Ti. When 
D=0.5/im snapshots are taken at time (b) SOT;, (c) 5000T;, 
(d)12000T,. 

can be further simplified when both Vg and Vp are very 
close to light speed, which means Wetc/t ^ {vp — Vg). The 
proton dephasing length is: 

Ldp^ iT^{aonc/ne)^/'^Xi. (4) 

In case of Li > j{aonc/neY^'^\i the laser pump deple- 
tion length is greater than the dephasing length , the 
maximum energy proton beam can be given as 

W,nax - ^{alnc/ne)m^c^. (5) 

It implies that the maximum proton energy in the snow- 
plow regime scales with Og and n~^. This is quite an 
efficient acceleration scaling law. 

The snapshots of electron density, proton density, and 
electrostatic field are plotted in Figl5] The proton beam 
has been continuously accelerated until laser pulse is 
completely depleted at t=12000Ti. When the proton 
beam approaches the laser front, the beam loading causes 
a substantial reduction of the electrostatic field as shown 
in Fig 121c) and^d). The theoretical values of the lon- 
gitudinal electrostatic field = 14.9 in the unit of i?o, 
the dephasing length Ipd = 11125/Ltm,the maximum en- 
ergy of proton beam Wmax = 532 GeV and the etching 
velocity Vetch = 0.0036c, respectively. They are consis- 
tent with the simulation results. Figl3] shows the energy 
spread (FWHM ) of the beam is less than 20 %. The 
maximum proton energy and averaged one are 540GeV 
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FIG. 3: (Color Online) Longitudinal phase space for protons 
Px/rriiC (a) and energy spectrum of trapped protons (b) from 
simulations at t=12000Ti. 

and 437GeV, respectively. There are 1.72% of the to- 
tal protons located inside this energy window. The pro- 
ton beam is compressed in the phase space and the final 
quasi-monoenergetic beam is shorter than 100 fj,m. 

As FigUl^a) shows, proton beams move backward with 
respect to the snowplough reference since the laser group 
velocity is greater than the proton beam. As a result, 
the distance between proton beam and laser front in- 
creases to reach the maximum li„j in the initial injection 
stage. Later it begins to decrease due to the erosion of 
the laser pulse front. The slopes of the solid lines repre- 
sent the etching velocities of the laser pulse for different 
gas densities. The maximum longitudinal electrostatic 
field, dephasing length, and maximum proton energy are 
plotted versus the density in Fig|4ljb,c,d). It shows the 
maximum energy of the proton beam is increasing with 
plasma density decreasing. The dotted lines indicate that 
the pulse depletion length is shorter than the dephasing 
length if the plasma density is smaller than O.Olric. In 
order to increase the proton energy, we may increase the 
intensity of lase pulse and reduce the gas density, which 
can increase both the dephasing length and pump deple- 
tion length. 

The density ranges over more than three orders of 
magnitude in ID simulations, it is difficult to do a full 
scale 2D simulations (e.g. TeV protons). A small sim- 
ulation box 40xlOOOA^ is chosen with a spatial resolu- 
tion of 5 cells/ A/ in y plane and 20 cells/ A/ in x plane. 
Each cell contains 4 particles for each species of gas 
plasma, 400 particles for foil plasma. The foil is the 
same as in ID simulations while the gas density is 0.2nc. 
The laser pulse is temporally and transversely super- 
Gaussian, / — /oexp(— (r/ro)^ — [{t — to)/T]^), where 
ro=10 /im, to=55 fs, and t=55 fs are taken. FigE] shows 
the snowplow layer, electrostatic field, and phase space are 
similar with ID results while the maximum proton en- 
ergy reaches 50GeV that agrees with ID simulations. 
Since protons gain energy mainly from the electrostatic 
field [IJ], the return current and magnetic field in 2D 
simulations are not important for proton acceleration. 
Furthermore, the lost fraction of electrons in snowplow 
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FIG. 4: (Color Online) (a)Distance between proton beam and 
the laser pulse front. It varies with time for different gas densi- 
ties, where blue, black, red and green curves represent the gas 
densities of O.Olric, 0.025nc, O.OSnc and O.lnc, respectively; 
(b)Longitudinal electrostatic field ei5a;/mea;iC;(c)Dephasing 
length; (d) Maximum proton energy. Here the solid lines 
are theoretical curves and simulation results are marked by 
different symbols(e.g. star, diamond and squre). 




FIG. 5: (Color Online)2D simulation results at t=990 
Ti. (a) Electron density distribution; (b)Electrostatic field on 
the axis ;(c)Proton phase space; (d)Proton spectrum. 

layer can be compensated from gas plasma and snowplow 
layer is still kept at the end of simulations. The disad- 
vantageous multi-dimensional effects (e.g. hole boring, 
Rayleigh- Taylor instability) in laser foil interactions ^] 
are not observed in our simulations. These effects lead to 
very stable proton acceleration in the multi-dimensional 
situation. 

In conclusion, snowplovif ion acceleration, using an 
ultra-relativistically intense laser pulse irradiating on the 
combination target, is presented. A ID model that can 



predict the proton dephasing length, pump depletion 
length, and maximum proton energy is derived, which 
is verified by ID and 2D PIC simulations. It shows 
that sub-TeV quasi-monoenergetic proton bunches can 
be generated by a centimeter-scale laser wakefield ac- 
celerator, excited by a laser pulse with the intensity of 
10^3 w/cm^ and duration of 116fs. The final proton 
beam is shorter than 100 /im and may be used to excite 
the plasma wakefield to accelerate electrons to hundreds 
GeV, as Cadwell proposed [l7| . 
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